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^ \ Abstract 

The anomalous couplings of the top quark and the Higgs boson have 
been studied in an effective theory deduced from the minimal supersym- 
, metric extension of the standard model (MSSM) as the heavy fields are 

| integrated out. Constraints on the parameters of the model from the 

^ ■ experimental data of Rb = T(Z — > bb)/T(Z — > hadrons) are obtained. 
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1 Introduction 



The standard model (SM) has been very successfull phenomenologically, but nevertheless it should be 
considered just as an effective theory valid for physics at the electroweak (EW) scale. In higher-energy 
regimes new physics beyond the SM must exist. Irrespective to what this new physics might be, it should 
be able to give a satisfactory answer to the most fundamental open question of the electroweak physics, that 
is, it must explain the origin of the electroweak gauge symmetry breaking [1, 2]. In the SM this is arranged 
through the spontaneous symmetry breaking mechanism by introducing a doublet of scalars with a nonzero 
vacuum expectation value (VEV). This mechanism, despite its simplisity and economy, has well known 
problems, which has enforced theorists and experimenters to look for new physics beyond the SM. Among 
the possible ways of extending the SM, supersymmetry is considered as a particularly attractive one. The 
minimal supersymmetry extension of the standard model (MSSM) provides an appealing solution to the 
gauge hierarchy problem by guaranteeing the perturbative stability of the theory from the electroweak 
scale to the Planck scale. 

The MSSM contains two complex Higgs doublets, denoted by H u , H d and assigned with opposite 
hypercharges Y B (H U ) = —Y B (H d ) = 1. There are altogether four neutral scalar degrees of freedom, three 
of which correspond to physical scalar fields. In the case where CP is conserved one can define two CP- 
even neutral Higgs fields, H, h, and one CP-odd neutral Higgs field, A. The present experimental bounds 
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on the Higgs boson masses set strong restrictions on the parameter space of the CP-conserving MSSM 
[3]. Radiative corrections to the lightest CP-even Higgs boson mass have been computed by using the 
renormalization group equation (RGE) method and diagram technique [4], and the resulting upper bound 
is 135 GeV, which is not much above the present experimental lower bound of 95GeV (95% CL). 

The possibility of CP violation makes the situation drastically different. There are three main sources 
of the CP violation in the MSSM Lagrangian. The first one is the well known \x parameter of the super- 
potential, which is in general complex. The second source is constituted by the soft mass terms of the 
SU(3) x SU{2) x U(l) gauginos. The third source are the phases of the soft supersymmetry-breaking 
mass terms of scalar fermions and of the soft trilinear couplings, which are presented by the matrices 
m Q u d l r an d \d ei respectively. Actually, only the off-diagonal elements of the soft mass matrices can 
be complex due to the hermiticity of these matrices. The matrices A [JC E in contrast can have complex 
phases also in their diagonal elements [5]. Not all the phases of these soft SUSY breaking parameters are 
physical and lead to the violation of CP parity. The physical CP phases are restricted by experimental 
observations, the most rigorous constraints originating from the measurements of the electron and neutron 
electric dipole moments (edm). The present upper limits for these edms are d e < 4.3 x lCT 27 e • cm [6] and 
d n < 6.5 x lCT 26 e • cm [7], respectively. Also the edm of P* 99 is quite accurately measured, the present 
upper limit being d 199 < 9. x 10~ 28 e • cm [8]. 

It has been demonstrated that the MSSM can be consistent with these constraints in some regions 
of the parameter space when suitable cancellations between different contributions occur [9] or when CP 
violation effects are associated with the third generation of squarks only [10]. The mixing of neutral Higgs 
bosons in the latter scenario is analyzed in [11, 12, 13, 14]. It is found that the CP-violating phases and 
large Yukawa couplings of the third generation fermions can lead to large mixings among the neutral Higgs 
bosons as a consequence of radiative effects. These mixings can drastically change the couplings between 
the neutral Higgs bosons and quarks and between the neutral Higgs bosons and gauge bosons, as well as 
the self-couplings of the Higgs fields. One consequence of this is that the experimental lower bound on the 
lightest neutral Higgs mass is relaxed to 60 GeV, while the predicted upper limit for the lightest Higgs 
boson mass remains about 135 GeV. 

If the new physics scale is much higher than the EW scale, one would have at the EW scale a great 
number of higher-dimensional operators Oi (dim{Oi) > 4) induced by the beyond-the-SM physics [15, 16, 
17, 18]. The resulting effective Lagrangian is of the general form 



Here Co is the SM Lagrangian, Cj are Wilson coefficients, and fi NP the energy scale of new physics. 
The Wilson coefficients are in general dependent on the new energy scale, but in addition to this all the 
higher-dimensional operators in C eff have a common suppression factor 1/ n 2 Np . 

In this paper we shall study anomalous couplings (to use the terminology of [19, 20]), i.e. the couplings 
not present in the SM Lagrangian Co, between the lightest neutral Higgs scalar (h) and the top quark 
induced by the new physics of MSSM. We assume that the other Higgs bosons, as well as all supersymmetric 
particles, are much heavier than the lightest neutral Higgs particle, so that the corresponding fields can be 
integrated out. A well known fact is that the masses of the other two neutral Higgs bosons are approximately 
equal to that of the charged Higgs boson (H + ) under the condition rn H+ ^> m h , and hence one can consider 
the lighter Higgs doublet as the SM Higgs field and integrate out the heavier Higgs doublet. 

Our presentation is organized as follows. In Section II, the notations adopted in this work are intro- 
duced. In Section III we shall describe the method of obtaining the Wilson coefficients by integrating out 
the heavy degrees of freedom in the full theory. The numerical analysis of the constraints on the parameter 
space from the present experiments, especially by the Rb data, is given in Section IV. Section V summarizes 
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our results. Some lengthy formulae, such as the expressions for the Wilson coefficients and the loop integral 
functions, are collected in Appendices. 



2 Preliminaries 

The most general gauge invariant superpotential, which retains all the conservation laws of the SM, is 
given by 

w = miKH + ^K'K 1 )^ + ^h J D I KQ I j t)J + ^K'KQW- ( 2 ) 

Here H u , H d are the two Higgs superfield doublets, Q 1 and L 1 are the doublets of quark and lepton 
superfields, and U , D 1 and R 1 are the singlet superfields of u- and d-type quarks and charged leptons, 
respectively (1=1, 2, 3 is generation index, i,j = 1,2 are SU(2) indices). Yukawa coupling constants are 
denoted by h L , h v D . The breaking of supersymmetry happens through the so-called soft terms, which are 
in the most general case given by 

C soft = -ml W*IP -ml H^H i -m 2 J ll*Li-m 2 TJ R I *R J -m 2 J1 Ql*Qi-m 2 J1 U I *U J 

-m 2 Dii b h b J + (uuXbXb + m 2 X i A X i A + m 3 X a G X% + h.c.) + [ £iJ m^ffff + eij A J J WL] R J 
^A^H\Q)D J + e^A^WQ 1 ^ + h.c] . (3) 

Here A% (a = 1, 2, ••• 8), (i = 1, 2, 3) and X B denote the 577(3), SU(2) and U(l) gauginos, 
respectively, and A l ^ D L are coupling constants of the unit of mass. 
Let us define scalar doublets $ and &h as follows: 

(:.)-(;?)(!)■ 

where H d = i^H* and c = cos/3 ,s = sin/3 with tan/3 = v u /v d , the ratio of the VEVs of H u , H d . 
With this definition $ is identified as the SM Higgs doublet, consisting of Goldstone bosons and a physical 
neutral Higgs field. More explicitly, one can write the two Higgs doublets as 

where G°, C + denote the Goldstone bosons, H® and H® are the neutral Higgs fields, H + and A are the 
physical charged Higgs and CP-odd neutral Higgs bosons, respectively, and v = ^jv\ + v 2 d = 246 GeV. 
At the electro-weak scale, the two physical CP-even neutral Higgs fields are obtained through the mixing 
between the fields Hf and H$. The masses of the physical Higgs bosons are given by 



m 2 = 

even 



(in the basis (H?, Hfj) T ), 
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S l3 C /3 
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m H+ = + • ( 6 ) 

t> l3 L l3 
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In the limit m H 3> the two doublets $ and decouple, the former remaining light and the latter 
being associated with a large mass m 2 TT = m 2 Ti ls R c„. 

H -"12 " ' 

In the following, we will use the four-component spinor representation for fermions. From the two- 
component Wyel spinors ip Qi , ip Ui , tp Di , ip Hu and ip Hd , we form the following four-component Dirac 
fermions: 

with tp Hd = {ia 2 ^jip Hd . Similarly, for the SU(3) x SU{2) x 17(1) gauginos Xq, X l A , \b can we define the 
following four-component Major ana spinors: 

Diagonalizing the soft mass terms is done with help of the sfermion mixing matrices Z QU D defined as 

Z l m l Z u = K > 



z K z o = K . (9) 

where the matrices rh^ pfl on the right-handed side are diagonal. 

Finally, we will benefit in our calculations from the following rearrangement identities of the SU(2) 
group indices: 

l aa > <8> 1/3/3' = ^ { 1 "/3 ® V"' + X] ® ^'a' } ' 

a 

a 

7^ tabctabd = 2b c d ■ (10) 
a, ft 



3 The Higgs-top anomalous couplings 

In this section we shall discuss the anomalous couplings of the top quark and Higgs bosons. Considering 
the suppression of the new physics energy scale, we just keep operators up to dimension-six in the effective 
Lagrangian Eq. (1). The top-Higgs anomalous couplings of interest can be classified into three types: the 
anomalous couplings involving a left-handed quark, the right-handed top quark and Higgs boson (O tq<s> ), 
the couplings between the Higgs boson and a left-handed quark (0 # ), and the couplings between the Higgs 
boson and the right-handed top quark (O ttt ). After the EW symmetry breaking, these operators produce 
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not only corrections to the effective couplings Wtb, Xti, Xbb (X = 7, Z, H), but also induce anomalous 
couplings such as jHti, ZHti. All those effects may be detectable at the Next Linear Collider (NLC) and 
at the Tevatron. 

In the following subsections we will give the explicit expressions for the contributions of supersymmetric 
particles and the heavy Higgs boson doublet to the effective operators mentioned above by deriving the 
relevant Wilson coefficients. We will give our results in terms of the following loop integral functions: 

(4tt) 2 f d A q (V 



^•'^ i J {2^Y^- Xa y^-x b f ' 
c < x x x) = ^l f^L 

V {x x x x) = ^lf^L W 

> MmK a ' 6 ' c ' d) i J (27r) 4 (q 2 -x a y(q 2 -x b ) k (q 2 -x c Y(q 2 -x d ) m ' 1 ' 

The explicit expressions of these are given in Ref. [21]. 

3.1 The anomalous couplings O tq<t 

This class of operators includes the CP-even operators 

O tq ^ = (^^)(q L t R ^ + m R q L ), 



O tq , 2 = q L (D,t R )D^ + (D^)\Dn R 



O tq * 3 = (D ll q L )(iyt R )* + &(D li t R )(D*q L ) , 
O tq * 4 = (D^q L )t R (D^) + [D^)h R (D»q L ) , 
O tq *s =i(D^q L )a^t R (D^) +i(p v ^ft R a^(p il q L ) , 



(12) 



where the covariant derivative is given by = <9 M — \g 3 T G ^ — ^g 2 aa W^ — ^g 1 Y B B ll . The CP-odd 
counterparts of these operators are 



O tq * 6 = (&$)(q L t R *-m R q L ) , 



O tq *r = Q L (D,t R )D^ - (D^)\DH R )q L , 
O tq * 8 = (D^ L )(DH R )$ - &(D,t R )(D 



£> tq * 10 = i(D^q L )a^t R (D^) - i(D u ®) k R a^ (d^) . (13) 
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For the CP-even operator O tq91 and for the corresponding CP-odd operator O t 96 , nonzero contribu- 
tions to the Wilson coefficients originate from the Feynman diagrams shown in Fig. 1, and they are given 
by 

<U = + g 2 2 )lm(h33) !^izjj , (14 ) 

where x ff = m2 H /^ 2 NP - I n the full theory, the Feynman diagrams that induce the anomalous couplings 
O tq<s>1 , O tq<Sf6 should also include diagrams involving virtual SM fields. However, these diagrams have no 
contribution to the Wilson coefficients after matching the effective Lagrangian Eq. (1) with the Lagrangian 
of the full theory (MSSM) (see below for more details). 

For the CP-even anomalous operators O tqii (i = 2, 3, 4, 5) and the CP-odd anomalous operators 
O tqi>i (i = 7, 8, 9, 10), the derivation of the Wilson coefficients leads to relatively tedious calculations. 
In Fig. 2 we show the Feynman diagrams, which induce nontrival contributions to the Wilson coefficients 
after matching the amplitude of the effective theory with that of the MSSM. In these diagrams, the black 
blobs represent the self-energy diagrams of q L q L , i R t R , and <J>t<J> (&h) (Fig- 4). 

The matching procedure, to which we refer above, is extensively applied in the derivation of the effective 
Lagrangian in the hadron physics, especially in the application of the effective Lagrangian to the rare B 
decay [22] . The main idea of this procedure is the following. We derive the amplitude corresponding to the 
relevant Feynman diagrams both in the full theory and in the effective theory. In both derivations we only 
keep the momenta pi of external particles to the second order. Through a comparison of the amplitudes 
of the full theory and the effective theory we then obtain the Wilson coefficients of interest. 

For a demonstration, let us consider the first diagram of Fig. 2. In the full theory we can write the 
amplitude corresponding to this diagram as 

A FT (p,q) = -— l—a c^fht h^hO ([A + l + ln^l(V$k„ V fl + In ^f) (q. $)q 2 t R 

Here A = i — j E + In 47r denotes the ultraviolet divergence (D = 4 — 2e is the time-space dimension in the 
dimensional regularization scheme), fi NP is the scale of new physics, and p and q denote the four-momenta 
of the external particles t R and <£, respectively. In the full theory the light fields and the heavy fields 
co-exist in the Lagrangian. When the heavy degrees of freedom are integrated out and the light fields are 
treated as massless, infrared divergences are encountered. They are regulated by the parameter m q . 

The amplitude of the corresponding Feynman diagram in the effective theory, presented in Fig. 3, is 
given by 



= -2uJ^^i( h MU A - 5 + ^ )M)' 2( - • (16> 



In the operators of the effective theory, only the light fields exist, and the Wilson coefficients do not 
depend on their masses. As in the full theory, an infrared divergence emerges here, and it is also regularized 
by the parameter m q . As expected, the infrared divergence appearing in the effective theory is the same 
as that appearing in the full-theory. By matching the amplitudes Eq. (16) and Eq. (15), one gets rid of 
the infrared divergence. After this matching step, we can present the amplitude in its final form: 
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(P, 9) = -jj^s^ 



(47T) 



(h f D h D h: 



:« 



A + 1 + In 



1 + ln^f 



'i:( 1+ln 



/v 

m 2 „ 



q L *)q-(p + q)t R + -^-(q L *)l>- (M f/)/ f 



(17) 



The first term in the parenthesis of Eq. (17) contributes to the renormalization of the Yukawa cou- 
plings h u , and it is irrelevant to our present discussion, taking into account the approximation level we 
work on. For those diagrams where the inner lines are supersymmetry particles, the Wilson coefficients 
of the anomalous couplings can be directly read from the amplitudes, because we integrate out all the 
supersymmetry fields in the effective theory. 

Now we will turn to show how to obtain the contributions of the self-energy diagrams to the anoma- 
lous couplings. As mentioned above, there are three possible self-energy diagrams that contribute to the 
coefficients indirectly, namely the self-energy corrections to the q L q L , i R t R , and Higgs doublet currents. 
For a fermion, the renormalized fields are defined by 



fO _ 72 f 

1 

fO _ 72 t 

J R,i «.'3 J Sj ' 



(18) 



where i,j are generation indices, /° ., f R . are the left- and right-handed bare fields, respectively, f Li , f Ri 
are the corresponding renormalized fields, and Z L R are the wave function renormalization constants. Ig- 
noring the fermion masses, we can write down the counter terms for the fermions in Eq. (18) as follows: 



^' c (p) = (Z%ZI M - S^j = \(5Zl i3 + 8Z LM )t , 



(19) 



where p denotes the external momentum of the fermion. In the full theory, we express the bare self-energy 
of the fermions as 
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ij 



5 +A L +B L p 2 U, 
S ij +A R + B R p 2 U 



(20) 



where the first term 5 tj represents the Born approximation part and A L,R , B L,R originate from raditive 
corrections. From Eq. (19) and Eq. (20) one finds the following form for the renormalized self-energies: 



6.. + -Uz\ .. +5Z r ..) +A L +B L p 2 U . 
5+1 (SZI +5Z R )+A R + B R p 2 U 



(21) 



The explicit expressions of the renormalization constant SZ L R depend upon the renormalization scheme, 
i.e., the renormalization conditions. Instead of the often-used renormalization schemes, i.e. the minimal 
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subtraction scheme (MS) or the modified minimal subtraction scheme (MS), we adopt here the physical 
renormalization conditions 





=0 


o, 


yRf 1 

ijJ R,j\f- 


=0 


o, 






/ L,i 


1 yRf 

ij J R>i 


1 j6=0 


~ f R,i 



(22) 

The first two conditions mean that the renormalized fields satisfy the equations of motion of free particles 
(for massless fermions this is a trivial constraint), and the last two conditions set the residue of the 
propagators at the pole equal to unity. In fact, this scheme is just the on-shell renormalization scheme 
often used when calculating radiative corrections to electroweak processes [23]. Of course, for high energy 
processes we can ignore the fermion mass in our approximation. Using the condition Eq. (22), we achieve 
the renormalized fermion self-energies: 



= B L P H , 

^ = ■ (23) 

We can attribute those terms to the contributions of the high dimension operators q L (jPj Q L , i R (jP^j t R ■ 
After the matching of the full and effective theories, there is no contribution to the operators of our 
interesting given in Eq. (12) and Eq. (13) from the fermion self-energy diagrams. 
For the Higgs sector, the bare self-energies are given as 

km) = d ** + i 1 + E **y + f **p 4 . 

K* H (P 2 ) = D ** H + E ** H P 2 + F ** H P 4 > ( 24 ) 

where p denotes the momentum of the external particle. In Eq. (24), D, E, F are standard integral 
functions that appear in radiative corrections. For the renormalization of the Higgs boson wave function 
and mass, we request the renormalized boson self-energy to satisfy the conditions 



£*JP 2 )I P2 . =0, 

1 

V 



-2^ $ (p 2 )| =0, 



5W/)| p2=0 =0, 

S„ AP% 2=m2 =0- (25) 



H 



It is easy to find the renormalized Higgs field self-energies which meet the conditions of Eq. (25): 

S„(P 2 ) = ^*P 4 , 

K, H (P 2 ) = F„ H P 2 (P 2 -0- (26) 
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The function is attributed to the contribution of the high dimensional operator ^{d^D^ <£. After 

the matching procedure, this piece will not contribute to the operators in (12) and (13) which we are 
interested in. In fact, after the matching the only nonvanishing contributions from the self-energy diagrams 
to these operators originate from the integral function , because we integrate the heavy Higgs doublet 
out in the effective theory 

After these preparations, we can now derive the Wilson coefficients of the operators O t iH (i = 
2, 3, 4, 5) and O tq9i (i = 7, 8, 9, 10). For clarity, we present their lengthy expressions in Appendix 
A. 

3.2 The anomalous couplings O t<s> 

This class of anomalous couplings includes the effective operators 

O tM = i(*t*) (i R Y(D^t R ) - {Dj- R )rt R ) , 

0„ 3 = i^D^ + (D^)i R ^t R , (27) 

where the operators O ti2 have even CP-parity and O t93 has odd CP-parity. In Fig. 5, we present 

those Feynman diagrams, which induce nontrivial contributions to the Wilson coefficients when matching 
the amplitude obtained in the effective theory with that in the full theory (MSSM). The ensuing Wilson 
coefficients are collected in Appendix B. 

In the full theory, we also include the 1PI diagrams depicted in Fig. 6, where the gray blobs represent 
the corresponding diagrams of Fig. 2. However, the contributions from these diagrams disappear as a 
result of the matching of the effective theory and full theory amplitudes. In order to demonstrate this, 
let us consider an example. From the subsection 3.1 , we find that the contributions of the subdiagram 
(framed by the dashed lines) in Fig. 7(a) induce the following term to the effective Lagrangian: 

1 5 

^ eff = 2U 2 ^ (Ptl** + ^t,*(5+a)) iPtqfc + ^tqt(5+ a )) ' ( 28 ) 



np a=2 



where 



+ 2 ^ A Z C rJ X ^ X ^ X Uj)] > 
^tq*3 = 247T 2 ^ 1 ^ %Vl S ' 3 ^'uj^ir.i( x ^^ x i' x u I ) ~ c /3 A U jCiy i ( x tJ.i x n x u I < 

C Ui = ~i§~2#lS K/^S^i'^'^j) - 2s ^Q^ x ^ x ^ x Uj)) + 2c p A ujQz( X ^ X i> X Q 1 ) 
C 't q *5 = ~^29i^2\ I Cl 12 (x fi ,x 1 ,x Ui ) , 
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C" 



t 9 *io 48^2 



(29) 



and x M = |/x| 2 /^ P , = m^//i^ p and x ; = |mj| 2 //i^ p (i = 1, 2, 3). The definition of the coupling 
constants and functions Qi(x,y,z) can be found in Appendix D. 
In the effective theory, the amplitude of Fig. 7(b) is written as 

+ (^ + CU)*-^9i-(91+Pi)*« 

where pi, gi, qi denote the four-momenta of the initial right-handed top quark and the Higgs bosons, 
respectively. In the full theory, the corresponding amplitude obtains the form 



96tt 2 h' 

"^121 ( X l* ' X l 1 X Uj ) + 2 ^"ll2 ( X M ' X l ' ) 

-tepQiix^x^x^in q! ■ (qi+Pi)t R 

2 ^i2i ( X it' X n x u I ) + ^ C f3 X v v /X M X i e ^ 1+Vm ^i3i(' C m' x i''%/) 



^ S B X ifi Y.i\ ( X M' ' X 



f 7 ' 



'*J*l + rfl 



Qi ■ Pit R 



+Ax 



Ul 



ft rp / rp 7y> ^ ( 1 ~ l~~ /-I ) Z*' / y rv» /-yj "\ 

/3 /av /a I L '113\ J 'A»' 1' E/// 



- C n 2 ( X m > x i > x c/ 7 ^ +^ W 1 >&]*«}> 



(31) 
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where (p^ and (fi (i=l,2,3) denote the CP-phases of the parameter [i and rrii, respectively. As already 
mentioned before, the first term of Eq. (31) is related to the Yukawa coupling renormalization in the full 
theory and it does not affect our computation. While matching Eq. (31) with Eq. (30), we find that the 
diagram does not contribute to the Wilson coefficients of the operators O t<s> . A similar conclusion is true 
also for the other 1PI diagrams in Fig. 6. 



3.3 The anomalous couplings O q<i 

This class of anomalous couplings includes the effective operators 



o q91 


= i(&D^-{D^*)q L >fq L , 




= i^a a D^ - (D^a a ^q L a a ^q L , 




= i(&*)(q L 'f(D ll q L )-(D ll q L )>fq L ) , 




= (q h a^{D^q L ) - (D^ L )a a ^q 




= i(^D^+(D^)^)q L ^q L , 


O 


= i(&a a D ix <!> + (D^<T a <S>)q L a a ^q L , 



(32) 

where the last two operators are CP-odd and the others are CP-even. The Feynman diagrams, which 
induce nontrivial contributions to the Wilson coefficients, are presented in Fig. 8. We collect the expressions 
for the Wilson coefficients of the corresponding operators in Appendix C. 



4 Experimental bounds on the Wilson coefficients 

At present, the most rigorous constraint on the Wilson coefficients considered in this work comes from 
the decay Z — ► bb. For an on-shell Z, one can write the general effective vertex Zbb as [19] 



r Zbb = -i- 



4s w c 



Vfa - A z w + ^S b z ( Pb - p- b )} , (33) 



where s w = sin(9 w , c w = cos# w , and p b , p b are the momenta of the outgoing quark and antiquark, 
respectively. For the operators listed in Eqs. (12), (13), (27) and (32), S z = 0. The vector and axial- 
vector couplings can be written as 

V z = V b Z '° + 6V z , 

A Z = A Z '° + SA Z , (34) 

where V Z '°, A Z, ° represent the SM couplings and SV Z , 5A Z are the new physics contributions. Ignoring 
the bottom quark mass, the lowest order theoretical prediction on the observable R b at the Z pole is given 
as 



_ T(Z -> bb) _ SM ( v z ^vZ + A Zfi 5A z SM , 

Rb ~ T(Z - hadrons) ~ ^ \ 1 + 2 ^^y7{^^ 
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With the Born approximation, we can obtain modifications to the couplings V b z , A b induced by the 
new physics operators #1 and O q<s>2 . Provided that there is no accidental cancellation between these 
contributions, the corrections are given by [20] 



0V b - 0A b - [C 9<E>1 + U q<Sf2 

' NP 



(36) 



where v denotes the VEV of the SM Higgs field doublet and m w is the lU-boson mass. From Eq. (35), 
we have 

" ~ " " (1 - R§ M )Rf M 2(V b z -" + A^) ' 1 ' 

The SM prediction on R b and the most recent experimental value are, respectively, given by [24]: 

Rf M = 0.21572 ± 0.00015, Rf v = 0.21664 ± 0.00065 . (38) 

If we attribute the difference of these two values to the new physics effects, we get a bound for the new 
physics corrections on the R b . At the la tolerance we obtain: 

0.00012 < AR b < 0.00172 . (39) 

Correspondingly, the bound for the Wilson coefficients is 

3.1 x 10- 4 < ^~C qH1+2) < 4.5 x 10- 3 (40) 

with C q9>(1+2) = C q<s>1 +C qi>2 . Using the same method, we can also analyze the forward-backward asymmetry, 
A h FB , of the decay Z — > bb. However, our theoretical result indicates that the present experiment data on 
this quantity set a weaker bound on the Wilson coefficients than R b . 

The other Wilson coefficients of the operators appearing in the Lagrangian are not constrained by R b on 
the Born approximation level. With higher-order approximations, those operators contribute to the gauge 
boson self-energies, and thus we can get for them only a rather loose bound with a significant uncertainty. 
We can also have loose bounds from the argument of partial wave unitarity [25]: 

-6-4 < C„, < 10.4 . (41) 

At present, there are no strong experimental constraints on the CP-odd couplings involving the top quark. 

It is well known that the MSSM contains in its general form unfortunately many 'new' free parameters 
in addition to the SM parameters. In order to simplify our discussion, we take the following assumption 
to restrict the MSSM parameter space: 

• All possible CP phases are taken to be zero or ir. A direct consequence of this choice is that there 
are no C-P-odd operators in the effective Lagrangian Eq. (1) 



All Yukawa couplings and the soft breaking parameters are flavor conserving, i.e., the mixing matrices 
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Under these assumptions, the parameters relevant to our discussion are the gauge coupling constants 
9ii 92 , 93, the higgsino and gaugino masses fi, mi, mi, m.3, the Yukawa couplings of the third gen- 
eration quarks and the corresponding soft breaking parameters h t = hP , hf, = hp , A t = A 33 , A^ = 
A 33 , and the square masses of the heavy Higgs boson doublet and the third generation squarks 
m 2 H , m 2 , m 2 , m 2 (I = 3). In our numerical analysis, we will disregard the loose bounds from partial 
wave unitarity on the Wilson coefficients C tqt>1 , C t#1 , C tqt>2 due to the large uncertainties mentioned above. 

Without losing generality, we assume m Q = m u = m D , A t = Af,, mi = 777,2 = 7773 in our numerical 
computations. Setting /j, NP = 1000 GeV, m H = 500 GeV, mi = 777,2 = 7773 = 500 GeV, A t = A^ = 
100 GeV, we obtain constraints set by Eq. (40) on the soft breaking parameters. In Fig. 9, we plot the 
values of m Q = m v = m D versus the parameter fi with (a) tan (3 = 2, and (b) tan (3 = 40, where the 
gray regions are allowed by the condition for v 2 C qiy(1+2) / fi 2 Np set by Eq. (40). From this plot we observe 
that the restriction set on the parameter space with tan (3 = 40 is more rigorous than that with tan (3 = 2. 
As tan (3 = 2, the contribution from the super symmetric box diagrams varies from negative to positive 
gradually, then tends to zero after its maximum as m Q = m u = m D increase from 200 GeV. When 
m Q = m u = m D > 1.3 TeV, the contribution is definitely less than 10~ 4 . Beside those box diagrams, C 91 
also receives a contribution from the heavy Higgs doublet. Under our choice about the parameter space, the 

Higgs contribution to the Wilson coefficient v 2 C q9>1 / '/z 2 p is proportional to 1/ tan 2 (3 — (m b /m^j tan 2 (3 . 
Taking the bottom quark mass m b = 4.5 GeV and top mass m t = 174 GeV, this contribution is about 
5 x 10 -4 for tan (3 = 2. As tan (3 increasing, the contribution of the heavy Higgs doublet is strongly 
suppressed and less than 10~ 7 for tan/3 = 40. The fact can help us understanding why very massive 
supersymmetry particles are allowed by the experimental bound for tan/3 = 2 (Fig. 9(a)), whereas the 
most part of the parameter space is excluded by the bound except a narrow band at the neighborhood of 
H = for tan/3 = 40 (Fig. 9(b)). 

In the figures discussed above, we have considered the la tolerance for the experimental data. Since 
the central value of the experimentally measured Rb is only about one standard deviation away from 
the SM prediction, this sets a lower bound on the C qiy{1+2) which is positive and very close to zero as 
shown in Eq.(40). Certainly, very massive supersymmetry particles are excluded by this condition in the 
large tan (3 case. In fact, considering the practical situation of the experiments, we may relax the lower 
bound on the C q<s>{1+2) to —5 x 10 -5 , while the upper bound remains unchanged (this is just only slightly 
beyond the standard deviation). In Fig. 10 we plot m Q = m u = m D versus ji by using the constraint 
-5 x 10~ 5 < v 2 C qiH1+2) /fi 2 Np < 4.5 x 10" 3 . One can see that the allowed parameter region is drastically 
enlarged in comparison with the case of the strict la tolerance for the large tan (3. 

Now, we discuss the operators #1 2 corrections to R b in the MSSM. Taking \i = rrii = 500 GeV (i = 
1, 2, 3), we plot AR b versus squark masses m Q = m u = m D with tan/3 = 2, 40 in Fig. 11. The gray 
band is the experimentally allowed region at the la tolerance. When the scalar quark mass is less than 
700 GeV, the supersymmetric box diagrams determine the leading contribution and results in a negative 
AR b , the corresponding parameter space is excluded by Eq. (39). As the parameters m Q = m u = m D 
increase, the supersymmetric contribution turns to be positive, then tends to zero after the maximum. 
With tan/3 = 2, the correction of the heavy Higgs doublet to the AR b is about 2.3 x 10~ 4 , and plays 
the leading role when m Q = m u = m D > 1.3 TeV. For tan/3 = 40, the total corrections from the Higgs 
and supersymmetric sectors to R b do not satisfy Eq. (39), because the contribution of the heavy Higgs 
doublet is strongly suppressed. For —fi = mi = 500 GeV (i = 1, 2, 3), the plot is similar to Fig. 11 
and not shown in the context. Taking m; = 500 GeV (i = 1, 2, 3), m Q = m u = m D = 500, 1000 GeV, 
and tan/3 = 2, 40, we present Ai? 6 versus the parameter fj, in Fig. 12. For m Q = m u = m D = 500 GeV 
(dot- and dot-dashed- lines), the corresponding parameter space is excluded by the condition Eq. (39) due 
to the negative supersymmetry contribution. With m Q = ra v = m D = 1 TeV and tan/3 = 2 (solid- line), 
AR b satisfies the condition Eq. (39) when \i > —700 GeV. As for the case m q = m u = m D = 1 TeV and 
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tan/3 = 40 (dashed-line), the new physics correction to AR b is excluded by the la tolerance experimental 
bound except the region neighboring /j = GeV. Finally, we investigate the new physics prediction on 
AR b with the assumption m Q = m u = m D = Choosing = 500 GeV (i = 1, 2, 3), tan/3 = 2, 40, 
we plot AR b versus the parameter m Q = m u = m D = |//| in Fig. 13. For the case tan/3 = 40, the 
correction to AR b exceeds the la tolerance experimental bound. As tan (3 = 2 and fi > 800 GeV, the 
new physics prediction on AR b satisfies this bound because of the relatively large contribution from the 
heavy Higgs. In those analyses the experimental bound with la standard deviation are adopted. After we 
relax the condition (Eq. 39) slightly, the more massive supersymmetry particles are also permitted by the 
corresponding experimental bound. 

Since the experimental data constrain the coefficients C qil 2 strongly, the operators O q<Sfl 2 have only 
negligible effects on the measurements at the proposed future colliders [19]. Other operators will produce 
the observable effects in the next generation colliders. In the associated production of the Higgs boson 
and top quark pair e + e~ — ► tt h, the CP-even operators will affect the energy and angular distributions of 
the final state particles [20]. Through the measurements of various distributions, such as da/dE t , da/dE h 
and da / d cos h (E t , E h denote the outgoing energy of the top quark and Higgs boson respectively, h is 
the angle of three-momentum of the outgoing Higgs boson with respect to the electron beam direction), 
we can obtain useful information about the operators. The constraints on the CP-odd operators can be 
obtained through measuring various CP violation observables in this process. In the process e + e~ — > tt, 
we can analyze the effects of the operators on various polarized top-quark production cross sections. On 
the other hand, more strict constraints on the supersymmetry parameter space will be set by more precise 
measurements on the widths of Z — > bb and the top quark decays. All of these will provide valuable 
information for the search of supersymmetry particles on the future colliders. 

It should be stressed that the above numerical analysis is performed under special assumptions about the 
MSSM parameter space. For example, we assume that all the parameters are real and flavor-conserving, the 
universal soft parameters are: m Q = m u = m D , A t = Ab, mi = mi = 771,3. I n a practical phenomenology 
analysis, those priori conditions should be dismissed. Nevertheless, it is quite obvious that the experimental 
data on Rf, set significant bounds on the parameter space even in a more general case than that we have 
considered here. 

5 The Summary 

We have considered in this work the anomalous couplings between top quark and Higgs boson induced 
by the MSSM when the heavy Higgs doublet and all supersymmetry fields are integrated out. An essential 
assumption made here is that there is only one neutral Higgs boson with the electroweak mass, the other 
Higgs particles are much heavier. We have derived the Wilson coefficients of the relevant higher dimensional 
operators in the ensuing effective theory. We have also studied numerically the constraints set by the 
experimental results for R^ = T(Z — > bb)/T(Z — > hadrons) on the parameters of the MSSM. 
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A The Wilson coefficients for the operators O t m (i = 2, • • • , 5, 7, • • • , 10) 
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i=l,2 



48vr 2 



^CCK^i'^) (42) 

with F 1 = 1, F 2 = 3, If = -T§ = I Y» = |, yj 3 = -| and A, = ^(a, - ^h,), A D = 
^(A c -^tan/3h D ), A E = ^ ( A E - // tan (3h E ) and A' f = A F + ^h F (F = U, D, E). In 
the above expression, the sum with the generation indices I, J is implied. 



B The Wilson coefficients for operators O tm (i = 1, 2, 3) 
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Just as in the appendix A, the sum with the generation indices /, J, K is implied. 

C The Wilson coefficients for operators O m (i = 1, • • • 6) 
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The coupling constants and loop functions 

The loop functions are denned as 
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Kj = V^(^[ Z Q( Z Q h V) 13 ] COS ^ + ^) " Im [ Z ?(4 ht J /3 ] sin ^ + ^)) . (* = !' or 2 ) 



A 



Af J = Re 



A 



Im 



Z ?(4 h t) ] "info*. + + Im^^ht 

^(4^Kh,) 3 3 + ^K h ^^)"4 J3 



/3 



cos(^ M + ipi) , (i = 1, or 2) 



A^ = Re((ht 4)3^f ) , A^ = lm((ht 4)3^^3) , 
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Re(z-(4ht )«) , A^=I m (^(4ht)-) ; 



c 2 Re 



s 2 Re 



(z u h u y 3I (z u h u h\ f zl) IJ (z u h 
(z D h D y 3 \z D h D tf D zl) IJ (z D h D ) 

{^) z \z D k D z Q ) 1J (z^l) 



J3 



J3 



cosipi-Im Zq 1 (z u A u Z q ^ Z j J" smipA , (i = 1, or 3) 



VM Re 



sin c^j + Im 

/J / . , x J3 



cos 



,73 



cos ipi) , (i = 1, or 3) 



I m 



.73 



^(R4K4) 3/ (^A D z Q ) /J (zt h t 



sin ^ 

73 



+Im 
Re 



h t D 4) 3/ (^A D Z Q ) /J (zt h t 



.73 



COS 



sin y? M 



+Im 



( h U Z Q ) 3/ (X K Z Q ) | J 4 3 ] [°f1 COS Vm + Sp V*~1 COS VI 



, ^{b* [(h^ Q ) 3/ (Z v k v Z Q ) \z { 



Im 



-Re 



cos ipn + Im 

S cos Vm + s /3 v^T cos ¥>i 
c /3\/^ sin VM - s^y^sin^i 



(h^ Q ) 3j (^A^ Q )^f 



sin 



(K z l) 3I ( z uK z Q ) u z q 3 } s ? (<? cos vm + ^ cos vi ) 



sin 



Re 



-Im 



(K Z l) 3I ( Z u K Z q) u Z l J3 ] s (sV^ sin VM - Sp^x-smifi) 
+Im[(h t y Zt)^(^A v Z Q ) /j Zt J3 ]s ;3 (c (3 ^:cos( / 9 M + s^^x7cos^) , (i = 1, 2) , 
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cos 



Im 



(ht^t) (z q k q z Q ) z l J3 



sin 



PM 



sin ^ 



r. 



r. 



+Im 

-Re 
+Im 

-Re 
-Im 



( h l Z l) 3I ( Z ^)u Z l J3 } COS ^ 
( h D Z i ) 3/ (X A d Z Q ) 7 J ^ J3 ] (s g jx; sin ^ - c (j y/x: sin ^) , 

C (i ( S /3 Sil1 Pm - C S V^i Sil1 V*) 



AT 



q,IJK 
B 

,IJK 



Xn 



Hef(h p Z Q ) s/ (Z,A,Z Q )t JJ (Z,A,Z Q ) J,f (Ztht) sf ) , (g = tf, Z> 



A" 



A: 



X ,Im((li ir Z Q ) s/ (Z,A,Z )t"(Z,A,Z g ) J,f (Zj fat )") , ( g = [/, D; X[/ = *J, X D = # , 

(*A*q)U*A*q) ^4*1 , (q = U, D; Xu = s 2 , X D = c 2 ) , 
^ lm [ Z l\ Z iK Z Q )\X Z ,K Z Q ) JK Z l K 3 ] , (9 = U, D; Xu = £ Xd = cj) , 



W^) 3 X Z A^)^A^L( Z ^) J ' (9 = ^ D; = ^ x ° = # • (46) 
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Figure 1: The Feynman diagrams inducing nontrivial contributions to the Wilson coefficients of the 



ators O. A1 and O. , fl 
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Figure 2: The Feynman diagrams inducing nontrivial contribution to the Wilson coefficients of the operators 
(i = 2, ■ ■ ■ , 5, 7 ■ ■ ■ , 10) in the full theory. 
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Figure 3: The Feynman diagram corresponds to the first diagram of Fig. 2 in the effective theory. 
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Figure 4: The Higgs self-energy diagrams which induce nonzero contributions to the Wilson coefficients of 
operators Eq. (12) and Eq. (13). 
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Figure 5: The Feynman diagrams which induce nontrivial contributions to the Wilson coefficients of the 
operators O uii (i = 1, 2, 3). 
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Figure 6: The one-particle- irreducible (1PI) Feynman diagrams which are related to the Wilson coefficients 
of the operators O t<Si (i = 1, 2, 3) in the full theory, the gray bulbs represent the diagrams of Fig. 2. 
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(a) 



(b) 



Figure 7: The 1PI Feynman diagrams related to the Wilson coefficients of the operators O tiSfi (i = 1, 2, 3) 
in (a) the full theory, (b) the effective theory. 
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Figure 8: The Feynman diagrams which induce nontrivial contributions to the Wilson coefficients of the 
operators O q9i (i = 1, • • • , 6). 
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Figure 9: The constraint from the anomalous coupling #1 2 set by the Rb experimental data with la 
tolerance, on the soft breaking parameters m Q = m v = m D versus the // parameter in the superpotential 
with mi = m 2 = m 3 = 500 GeV, A t = A b = 100 GeV, m H = 500 GeV and (a) tan (3 = 2; (b) tan (3 = 40. 
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Figure 10: Relaxing the lower bound to — 5x and keeping the upper bound unchanged as in Eq.(40), the 
constraint from the anomalous coupling O qi>1 2 on the soft breaking parameters m Q = m u = m D versus the 
/i parameter in the superpotential with m\ = m2 = = 500 GeV, A t = = 100 GeV, m H = 500 GeV 
and (a) tan (3 = 2; (b) tan (3 = 40. 
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Figure 11: Taking p = = 500 GeV (i = 1, 2, 3), A t = A b = 100 GeV, m H = 500 GeV, AR b versus 
squark masses m Q = m v = m D with tan (3 = 2 (solid line) or tan (3 = 40 (dashed line). 
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Figure 12: AR b versus the parameter [i with rrii = 500 GeV (i = 1, 2, 3), At = Af, = 100 GeV, m H = 
500 GeV, and (a) solid-line: m Q = m u = m D = 1 TeV, tan (3 = 2, (b) dashed-line: m Q = m u = 
m D = 1 TeV, tan/3 = 40, (c) dot-line: m Q = m u = m D = 500 GeV, tan/3 = 2, (d) dot-dashed-line: 
m Q = m v = m D = 500 GeV, tan (3 = 40. 
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Figure 13: Under the assumption m Q = m v = m D = \fi\, AR b versus the parameter jjl. The other 
parameters are taken as rrii = 500 GeV (i = 1, 2, 3), A t = Af, = 100 GeV, m H = 500 GeV and tan/3 = 2 
(solid- line) tan/? = 40 (dashed- line). 
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